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Abstract 

In this article we will present pure three dimensional analytic solutions for the Navier-Stokes 
and the continuity equations in Cartesian coordinates. The key idea is the three-dimensional 
generalization of the well-known self-similar Ansatz of Barenblatt. A geometrical interpretation 
of the Ansatz is given also. The results are the Kummer functions or strongly related. Our final 
formula is compared with other results obtained from group theoretical approaches. 
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To describe the dynamics of viscous incompressible fluids the Navier-Stokes (NS) partial 
differential equation (PDE) together with the continuity equation have to be investigated. 
In Cartesian coordinates and Eulerian description these equations have the following form: 

Vv = 0, 

Vi + (vV)v=i/Av-y + o (f) 

where v,p,p,is,a denote respectively the three-dimensional velocity field, density, pres- 
sure, kinematic viscosity and an external force (like gravitation) of the investigated fluid. 
(To avoid further misunderstanding we use a for external field instead of the letter g which 
is reserved for a self-similar solution.) In the following is, a are parameters of the flow. 
For a better transparency in the following we use the coordinate notation v(x,y, z,t) = 
u(x, y, z, t),v(x, y, z, t),w(x, y, z, t) and for the scalar pressure variable p(x, y, z, t) 



U X + Vy + W Z = 

I \ Px 

U t + UU X + VUy + WU Z = V[U XX + Uyy + U ZZ ) ~ ~ 

v t + uv x + vvy + wv z = is(v xx + v yy + v zz )-^ 

Pz 

w t + uw x + vw y + ww z = is(w xx + Wyy + W zz ) — + a. (2) 

The subscripts mean partial derivations. According to our best knowledge there are no 
analytic solution for the most general three dimensional case. However, there are various 

n 

examination techniques available in the literature. Manwai [1| studied the N-dimensional 
(N > 1) radial Navier-Stokes equation with different kind of viscosity and pressure depen- 
dences and presented analytical blow up solutions. His works are still 1+1 dimensional (one 
spatial and one time dimension) investigations. Another well established and popular in- 
vestigation method is based on Lie algebra there are numerous studies available. Some of 
them are even for the three dimensional case, for more see J2|. Unfortunately, no explicit 
solutions are shown and analyzed there. Fushchich et al. 3] construct a complete set of 
G(l, 3)-inequivalent Ansatze of codimension 1 for the NS system, they present 19 different 
analytical solutions for one or two space dimensions. They last solution is very closed to 
our one but not identical, we will come back to these results later. Further two and three 



dimensional studies based on group analytical method were presented by Grassi ^]. They 
also present solutions which look almost the same as ours, but they consider only 2 space 
dimensions. We will compare these results to our one at the end of the paper. 



I. E 



Recently, Hu et al. [5] presents a study where symmetry reductions and exact solutions 



of the (2+l)-dimensional NS were presented. Aristov and Polyanin jg] use various methods 
like Crocco transformation, generalized separation of variables or the method of functional 
separation of variables for the NS and present large number of new classes of exact solutions. 
Sedov in his classical work [\* presented analytic solutions for the tree dimensional spherical 
NS equation where all three velocity components and the pressure have polar angle depen- 
dence (9) only. Even this kind of restricted symmetry led to a non-linear coupled ordinary 
differential equation system with has a very rich mathematical structure. 

Beyond the NS system there are other important and popular PDEs which attract 
much interest and investigation. The applied methods are the same there, too. With- 
out completeness we mention some examples. For one dimensional cubic-quintic nonlinear 
Schrodinger equation a quite general self-similar type of solution ip(z, t) = u(z, t)exp[iv(z, £)] 
was applied where u and v are real functions [8j. The results are analytic solutions 
for an external potential with variable coefficients. A more general type of this Ansatz 
u(z,t) = A(z)U[T(z,t)exp(i<p(z,t))] was used with success to get chirped and chirp-free 
self-similar cnoidal solitary wave solutions [9[ for the same equation. Such solutions can 
be generalized for multi dimensional spatial coordinates. There are analytic solitary wave 
solutions available for the (3+1) dimensional Gross-Pitaevskii equation with the following 
Ansatz ip = u(x, y, z, t)R(t)exp[ib(t)(x 2 + y 2 + z 2 )} [lo| . 



From basic textbooks the form of the one-dimensional self-similar Ansatz is well-known 

Mm 

T(x,t) = ry^) :=r/W, (3) 

where T(x, t) can be an arbitrary variable of a PDE and t means time and x means spatial 
dependence. The similarity exponents a and /3 are of primary physical importance since 
a represents the rate of decay of the magnitude T(x,t), while ft is the rate of spread (or 
contraction if /3 < ) of the space distribution as time goes on. The most powerful result of 
this Ansatz is the fundamental or Gaussian solution of the Fourier heat conduction equation 
(or for Fick's diffusion equation) with a. — (3 — 1/2. These solutions are visualized on figure 
1. for time-points t\ < t%. In the pioneering work of Leray [3] in 1934 at the end of the 
manuscript he asks whether it is possible to construct self-similar solutions to the NS system 
in R 3 in the form of p(x,t) = ^P(x/VT^t) and v(x,t) = j=V ( x I 'VT^) ■ In 2001 
Miller et al. 1^] proof the nonexistence of singular pseudo-self-similar solutions of the NS 
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FIG. 1: A self-similar solution of Eq. Q for t\ < t-i. The presented curves are Gaussians for 
regular heat conduction. 

system with such kind of solutions. Unfortunately, there is no direct analytic calculation 
with the 3 dimensional self-similar generalization of this Ansatz in the literature. We will 
show later on that in our case the time dependence has the same exponents as showed above. 
Applicability of this Ansatz is quite wide and comes up in various transport systems 



12 



15 



171 ] . This Ansatz can be generalized for two or three dimensions in various 



ways one is the following 

»,„,*) = r«/ ) := r-/ (£±»±£) := r«/M (4) 

where F(x, y, z) can be understood as an implicit parameterization of a two dimensional 
surface. If the function F(x, y,z)— x + y + z — which is presented on figure 2. then it 
is an implicit form of a plane in three dimensions. At this point we can give a geometrical 
interpretation of the Ansatz. Note that the dimension of F(x, y, z) still have to be a spatial 
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FIG. 2: The graph of the x + y + z = plane. 



coordinate. With this Ansatz we consider all the x coordinate of the velocity field v x = u 
where the sum of the spatial coordinates are on a plane on the same footing. We are 
not considering all the R s velocity field but a plane of the \ x coordinates as an independent 
variable. The Navier-Stokes equation - which is responsible for the dynamics - maps this kind 
of velocities which are on a surface to another geometry. In this sense we can investigate 
the dynamical properties of the NS equation truly. In principle there are more possible 
generalization of the Ansatz available. One is the following: 



which can be interpreted as an Euclidean vector norm or L 2 norm. Now we contract all 
the x coordinate of the velocity field u (which are on a surface of a sphere with radius 
a) to a simple spatial coordinate. Unfortunately, if we consider the first and second spatial 
derivatives and plug them into the Navier-Stokes equation we cannot get a pure r\ dependent 
ordinary differential equation(ODE) system some explicit x, y, z or t dependence tenaciously 
remain. For a telegraph-type heat conduction equation both these Ansatzes are useful to 
get solutions for the two dimensional case jl7 ]. 

Now we concentrate on the first Ansatz (4) and search the solution of the Navier-Stokes 
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PDE system in the following form: 

u{x,y,z,t) = t } I p I, v(x,y,z,t) = t >g I ^ I, 

y, z, t) = t~ e h (^ X+ ^ c + , p(x, y, z, t) = m + ^ + ^ . (6) 

Where all the exponents a, (3, 7, 5, e, £, 77, 6* are real numbers. (Solutions with integer expo- 
nents are called self-similar solutions of the first kind, non-integer exponents mean self-similar 
solutions of the second kind.) The functions /, g, h, I are arbitrary and will be evaluated later 
on. According to Eq. ([2]) we need to calculate all the first time derivatives of the velocity 
field, all the first and second spatial derivatives of the velocity fields and the first spatial 
derivatives of the pressure. All these derivatives are not presented in details. Note that both 
Eq. (T5]) and Eq. ([8]) have a large degree of exchange symmetry in the coordinates x, y and z. 
Later we want to have an ODE system for all the four functions f(u)),g(u)), h(cu), l(u) which 
all have to have the same argument u. This dictates the constraint that (3 = 5 = ( = 9 
have to be the same real number which reduces the number of free parameters, (let's use 
the (3 from now onu= x +v+ z ) m From this constrain follows that e.q. u x = ~ v y = 
where prime means derivation with respect to u. This example shows the hidden symmetry 
of this construction which may helps us. For the better understanding we present the second 
equation of (T5]) after the substitution of the Ansatz ([6]). 

- at-^fico) - pr^n^u; + r 2a -Pf(u)f'(cu) + t-^-^M/V) + 

r e - a -Ph{uj)f{uj) = u3r a -^f"(u) - L!!_JM. (7) 

To have an ODE which only depends on u (which is now the new variable instead of time 
t and the radial components) all the time dependences e.g. have to be zero OR all the 

exponents have to be the same. After some algebra it comes out that all the six exponents 
a — ( included for the velocity filed (the first three functions in Eq. ()8])) have to be +1/2. 
The only exception is the term with the gradient of the pressure. There 77 = 1 and 9 = 1/2 
have to be. Now in Eq. Qj each term is multiplied by t~ 3 ' 2 . Self-similar exponents with the 
value of +1/2 are well-known from the regular Fourier heat conduction (or for the Fick's 
diffusion) equation and gives back the fundamental solution which is the usual Gaussian 
function. For pressure the rj = 1 exponent means, a two times quicker decay rate of the 
magnitude than for the velocity field. 
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FIG. 3: The KummerM 



I' l> ^llir ) function for c = 1, and f = 0.1. 
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FIG. 4: The KummerU 
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function for c = 1, and v = 0.1. 



Now we may write down the concrete form of the Ansatz ([6]) 

i(x, y, z, f) = r 1 / 2 / ( g+ t f /2 +Z ) = r^/M, y, *, f) = t-Wgiu,), 

w(x, y, z, t) = r 1/2 h(u), p(x, y, z, t) = t~H(u), 



(8) 



and the corresponding coupled ODE system 
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FIG. 5: The implicit plot of the self-similar solution Eq. (|17p . Only the KummerU function is 
presented for t = 1, c\ = 1, c% = 0, a = 0, c = 1, and v = 0.1. 
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/V) + <?V) + /*V) = o 

- + [/(a;) + <?(u,) + = 3i//» - 

--h(cj) - iw/i^w) + [/(w) + g(cj) + h(u)}h\u) = 3uh"(u) - + a. (9) 

Zi Zi j) 

From the first (continuity) equation we automatically get 

f(u) + g(u) + h(oo) = c, and f"(u) + g"(u) + h"(co) = (10) 

where c is proportional with the constant mass flow rate. Implicitly, larger c means larger 
velocities. From the second equation we can express — - and can substitute it into the third 
and fourth equation. After some algebra we arrive at 

+[/M + g (u>) + /iMKsV) - /V)) = o 



2 2 



(11) 
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Now inserting f"(u) = -g"(u) - h"(cu), f'(u) = -g'(cu) - h'(cu) and f(u) =c-g{u) - h(u) 
we get the final equation 



9i//» - 3(w + c)f'(u) + 5/( w ) -£ + fl = 0. 



(12) 



The solutions are the Kummer functions 

1 1 (oo + c) 2 



18| 



f(oj) = C\ ■ KummerU 



,,, 11 (w + c) 2 \ c 2a 

1-2- 0, /' +C2 '^ mmerM i-4'2'-6^J + 3-y 

(13) 



where Ci and C2 are integration constants. The KummerM function is defined by the following 
series 

Ma, 6,^ = 1 + — + y -£- + ... + ^-r 14 
6 (6) 2 2! (6) n n! 

where (a) n is the Pochhammer symbol 



(a) n = a(a + l)(a + 2). ..(a + n — 1), (a) = 1. 



(15) 



The KummerU function is defined from the KummerM function via the following form 



U (a, b, z) 



71 



sin{7ib) 



M(a,b,z) 



r(i + o-6)r(6) 



— z 



!_ 6 M(1 + a - 6,2 - b, z) 



T(a)r(2-6) 



(16) 



where T() is the Gamma function. Exhausted mathematical properties of the Kummer 



function can be found in 



1* 



Note, that the solution depends only on two parameters where the v is the viscosity, and 
c is proportional with the mass flow rate. Figure 3 and figure 4 show the KummerM and 
KummerU function for c = 1 and v = 0.1, respectively. For stability analysis we note that 
the power series which is applied to calculate the Kummer function has a pure convergence 
and a 30 digit accuracy was needed to plot the KummerU function, otherwise spurious 
oscillations occurred on the figure. Note, that for u = 6.5 the KummerM goes to infinity, 
and bj — > oo KummerU function goes to oo which is physically hard to understand, which 
means that the velocity field goes to infinity as well. 

The complete self-similar solution of the x coordinate of the velocity reads 

-1 1 ((x + y + z)/t l l 2 + c) 2 * 



u(x,y,z,t) 



t ^ 2 f{u 
t - 1 ' 2 



t 



-1/2 



Ci • KummerU 



4 '2' 



6/v 



r , • l\ iiiiniicrM j , -, — — - — — ~*~ - 

1 4' 2' 6u 



c 2a 

3 ~ y 



+ 



(17) 



9 



On figure 5 an implicit plot of Eq. ( ITT]) is visualized. The KummerU function was presented 
only, the used parameters are the following ci = l,C2 = 0,t = l,c = l,z/ = 0.1,a = 0. Note, 
that the initial flat surface of figure 2 is mapped into a complicated topological surface via 
the NS dynamical equation. The following phenomena happened, an implicit function is 
presented, we already mentioned that all the x + y + z = points considered to be the 
same. Therefore we got a multi-valued surface because for a fixed x numerical value various 
y+z combinations give the same argument inside the Kummer function. Unfortunately, this 
effect is hard to visualize. This can be understand as a kind of fingerprint of a turbulence- 
like phenomena which is still remained in the equation. An initial simple single- valued plane 
surface is mapped into a very complicated multivalued surface. Note, that for a larger value 
(now we presented KummerU() = 2 case) or for larger flow rate (c=l) the surface got even 
more structure. Therefore figure 5 presents only a principle. At this point we can also give 
statements about the stability of this solution, the solution the Kummer functions are fine, 
but for larger flow values a more precise and precise calculation of the solution surface is 
needed which means larger computational effort which is well known from the application 
of the NS equation. 

From the integrated continuity equation (/ = c — g — f) we automatically get an implicit 
formula for the other two velocity components 



v(x,y,z,t)+ w (x,y,z,t) = -t 1/2 
t - 1 ' 2 



c\ ■ KummerU I — , -, — — - — — — — — 

4 2 6u 

9n 1 

+ c(18) 



11 ((x + y + z)/^' 2 + cf\ c 2a 



For explicit formulas of the remaining two velocity components the two ODEs of ( ITTi) have 
to be integrated. For v(x,y,z,t) = t~ l l 2 g(uj) the ODE is the following 

- 3vg"(u) + g'{u) (-| + c) - ^ + F(f"(u), f'(co), f(co)) = (19) 

where F(f"(u), f'{uj), f{uj)) contains the combination of the first and second derivatives of 
the Kummer functions. This is a second order linear ODE and the solution can be obtained 
with the following general quadrature 



-d + / F(/», f(u), f(u))du ■ exp (^g^) 
' 2 / { ^ ~}du 



exp 



-oo 2 /4: + COJ 
3^ 



(20) 
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For the sake of simplicity we present the formulas of the first and second derivatives of the 

KummerU functions only 

d T ^ TT/ , . (uj + a — b) KummerU ( a, b,u) — KummerUia — 1,6, a;) 

— KummerU (a, o,ui) = (21) 

dui to 



and 

d 2 



-KummerU {a, b,u) = — 
duj 2 to 2 



a {uaU(a + 1, 6, u) — loU (a + 1, 6, u)b + loU (a + 1, 6, u) 



-aU(l + a, 6, u)b + U(a, 6, w)6 + U(a + 1, 6, w)6 2 - U(a + 1, 6, u)b}] .(22) 

Unfortunately, we could not find any closed form for v(x,y,z,t) and for w(x,y, z,t). Only 
v the x coordinate of the velocity v field can be evaluated in a closed form. 

As we mentioned at the beginning there are analytic solutions available in the literature 
which are very similar to our one. Fushchich et al. [3J present 19 different solutions for the 
full three dimensional NS and continuity equation. (For a better understanding we used the 
same notation here as well.) For the last (19th) solution they apply the following Ansatz of 



u(z,t) = v(y,z) 



?M + y w(zt) = } M v(tz) = l M 

Vt + V W[Z ' t} Vi ' K ' ] y/t 



(23) 



where uj = z/yt is the invariant variable. The obtained ODE is very similar to ours (j9J) 

h\u) + 1 = 
- l -(f(uj)+uf{uj)) + h{uj)f{u) = /», 

-(g(uj)+ug'(u)) + h(uj)g'(uj) = g"(uj), 

--(h(u)+uh'(u))+h(u)h'{u) + l'(u) = f"(u). 

/Li 



(24) 



The solutions are 

m 

h{u) 



-uj — c) 1 ^ 2 exp 



,3 



c) 1 ^ 2 exp 



13 

-6 ( 2 W - C) 
13 

■6 ( 2 W " C) 



1113 
W -12'4'3 ( 2 W - C) 

5 1 1,3 
" -12'4'3 ( 2 W - C) 



-uj + c 



-CUJ — UJ + C\ 



(25) 



where w is the Whitakker function, c and C\ are integrational constants. Note that the 



Whitakker and the Kummer functions are strongly related to each other [181]* 
w(k, h, z) = e~ 1/2z z 1/2+ ^KummerM(l/2 + y, - k, 1 + 2fi, z). 



(26) 
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More details can be found in the original work j^]. 

As a second comparison we show the results of [4J. They also have a modified form of 
(T2J) which is the following 

U u + cU x + U 2 U ly + U 3 U lz - u(U lyy + U U z) = 0, 

U 2t + U 2 U 2y + U 3 U 2z + 7Ty - v(U 2yy + U 2zz ) = 0, 

U 3t + U 2 U 3y + U 3 U 3z + tt z - v{U 3yy + U 3zz ) = 0, 

U 2y + U 3z + c = (27) 

where Ui, i = 1..3 are the velocity components Ui(y, z, t) and 7r is the pressure, c stands for 
constants, v is viscosity and additional subscripts mean derivations. After some transfor- 
mation they get a linear PDA as follows 

U u + k lV U ly + {a - k x z)U lz - v{U Xyy + U lzz ) = (28) 

it is convenient to look the solution in the form of 

U x = Y(y)T{z)$(t). (29) 

Note, that they also consider the full 3 dimensional problem, but the velocity filed has a 
restricted two dimensional (y,z) coordinate dependence. There are additional conditions but 
the general solution can be presented 

$ = ciexp(c 2 )t 

1 V 2 \ ,,fl 3 y 2 



Y = c 3 M|-c 4 ,-,|-j +!/ c 5 M^-c 1 ,-,|- 
T ^M[c^£\ + ,m(\-cJ-,£\ (30) 



where M is the Kummer M function as was presented below. The exact solution in {4] 
(4.10a-4.10c)* contains more constants as presented here. It is not our goal to reproduce 



the full calculation of j4| (which is not our work) we just want to give a guideline to their 
solution vigorously emphasising that our solution is very similar to the presented one. 
Note that in both results the arguments of the Kummer M function ( 1T3"|) and fl3U|) are 
proportional to the square of the radial component divided by the viscosity, additionally 
one of the parameters is 1/2. As a last word we just would like to say, (as this example 
clearly shows) that the Lie algebra method is not the exhaustive method to find all the 
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possible solutions of a PDA. 

In summary: We introduced and gave a geometrical interpretation of a three-dimensional 
self-similar Ansatz. We applied it to the three-dimensional Navier-Stokes equation in Carte- 
sian coordinates. The question of another Ansatze was mentioned briefly as well. Some part 
of the results could be written as Kummer functions. Unfortunately, some other parts of 
the results could not be written in closed forms. Further work is in progress, (we still have 
some hope) to learn something new from Eq. ( TT9|) . We compared our results with other 
analytic solutions obtained from various Lie algebra studies. The structure of the result - 
the implicit coordinate dependence of the Kummer function - was analyzed as well. We 
hope that even this moderate result can give any simulating impetus to the investigation of 
the Navier-Stokes equation. Our solution can have some real interest and can be used as 
a test case for various numerical methods or commercial computer packages like Fluent or 
CFX. 

The paper is dedicated to my first mathematics teacher " Sanyi Bacsi" . 
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